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INTRODUCTION 


In analyzing heating on vehicles such as the space shuttle orbiter or 
other advanced entry configurations , it is often necessary to calculate the 
heating rate at the stagnation point of these vehicles (which are usually 
three-dimensional) or to calculate a reference heating rate which can be used 
to nondlmensionalize the other surface heating rates experienced by the 
vehicle. The reference value usually chosen is the stagnation point of a 
sphere or other simple shape. 

These stagnation point heating rates can be obtained from a solution of 
the full Navler-Stokes equations or from a coupled solution of the Euler and 
boundary-layer equations. However, either of these procedures can become very 
expensive if a large number of calculations are required - such as over an 
entire entry trajectory. 

Reshotko (ref. 1) has developed a simple method of calculating the 
heating rate at general three-dimensional stagnation points which is 
applicable to cold wall flight conditions. Using this method, the heating at 
a general three-dimensional stagnation point can be related to the heating at 
an axisymmetric stagnation point by a simple algebraic expression which can be 
evaluated from the geometry of the stagnation region. 

Although several simple methods exist for calculating the heating rates 
at axisymmetric stagnation points, a method developed by Fay and Riddell 
(ref. 2) has been widely used and has been shown to yield reasonably accurate 
results compared with both experimental data and more detailed calculations 
(refs. 1 and 3). Combining the equation developed by Fay and Riddell for 
axisymmetric stagnation points with the equation developed by Reshotko 
discussed previously, a single simple algebraic equation is obtained for 
calculating the heating rate at general three-dimensional stagnation points 
which should be sufficiently accurate for most engineering calculations. 

The purpose of the present paper is to give a detailed description of the 
application of the method outlined previously for making stagnation point 
heating rate calculations for a vehicle during atmospheric entry. 




SYMBOLS 


stagnation point velocity gradients, 1/sec 

static enthalpy, ft^/sec^ 
dissociation enthalpy, ft^/sec^ 
total enthalpy, ft^/sec^ 
parameter defined by eq. (5) 

Lewis number 
mass fraction 
Mach number 

static pressure, Lb/ft^ 

Prandtl number 

heating rate, BTU/ft^-sec 

radius of curvature at stagnation point in x-direction 
(see Fig. 1), ft 

radius of curvature at stagnation point in z-direction 
(see fig. 1), ft 

velocity in x-dlrection (see fig. 1), ft/sec 
velocity along stagnation streamline, ft/sec 
velocity in z-direction (see fig. 1), ft/sec 
mole fraction 

compressibility parameter defined by eq. (21) 
parameters defined by eqs. (22) and (23), respectively 
ratio of specific heats for perfect gas 
viscosity, slug/ft-sec 
density, slug/ft^ 



Subscripts 


APPROX 

AXI 

BL 

e 

s 

sw 


w 


3D 


approximate solution 
axisymmetric 
boundary-layer solution 
boundary-layer edge 
stagnation point 
shock wake 
wall 

three dimensional 
free stream 


METHOD OF SOLUTION 


General Three-Dimensional Stagnation Points 


Consider the flow in the vicinity of a general three-dimensional 
stagnation point - that is, a stagnation point with different principal radii 
of curvature as shown in figure 1. Using the results from Reshotko (ref. 1), 
it can be shown that for a cold wall (which is usually the case in flight) the 
heating rate at a general three-dimensional stagnation point can be related to 
the heating rate at an axisymmetric stagnation point by the following simple 
expression 


^‘^s^ 3 D 1/2 ^‘^s^AXI 


( 1 ) 


where k is a parameter that is related to the two principal velocity gradients 
by the following expression (see fig. 1) and (qs)AXI based upon the 
small radius of curvature. 

^ (dwe/dz)s 

(dUg/dx)s 

Now if the x coordinate is always identified with the smaller radius of 
curvature Rj^, then the parameter k will always fall within the range 

0 < k < 1 (3) 
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dWo 

For two-dimensional flow R 2 «, ^ 

dz 

k = 0 (2D flow) 

dw dll 

while for axisymmetric flow R 2 = R^f f. = and 

dz dx 


k = 1 


(axisymmetric flow) 


Now using a modified Newtonian flow model, the following expressions for the 
two principal gradients can be obtained. 







(4a) 


(4b) 


Combining equations (3) and (4), the following is obtained 
k = Rx/Rz 


(5) 


Both Rjj and R 2 can be determined once the geometry in the stagnation 
region is known. This is discussed in detail in reference 14. 


Fay and Riddell (ref. 2) have developed the following expression for the 
heating rate at an axisymmetric stagnation point for a fluid in thermodynamic 
equilibrium 


^‘^s ^AXI 


0.76 

Tp?y° 




X 



( 6 ) 


This expression has been shown to yield results that are in good agreement 
with both experimental data and more detailed boundary layer calculations 
(refs. 2 and 3). 
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Combining equations (1) and (6), the following expression is obtained for 
the heating at a general three-dimensional stagnation point 


) = 0.76 

^s-'axI 


0. 6 






(7) 


where the velocity gradient (du^/dx)o is given by equation (Aa). 

( 1 + k 
— 

2 

in equation 7 is to yield an effective velocity gradient whict 
of the two principal velocity gradients. 



Stagnation Point Properties 

The flow approaching a stagnation point first passes through a normal 


shock wave and then compresses until 

The pressure (Pg^) and enthalpy 
shock wave can be expressed in terms 
(ref. 5) 



the velocity is zero (see fig. 2). 

(hsw) on the downstream side of the 
of the density (pg^) as follows 

( 8 ) 

1 (9) 


For an ideal gas, the density (pgw) on the downstream side of the shock is 
given by (ref. 6): 


Pgw = 

Poo (Y - mi + 2 


( 10 ) 
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Using this result equations (8) and (9) can be solved directly for the 
pressure and enthalpy. 

For equilibrium air chemistry, no simple explicit expression (such as 
eq. (10)) can be obtained for the density. Thus equations (8) and (9) must be 
solved Iteratively with an equilibrium equation of state in the functional 
form 


P = p(P»h) (11) 

to obtain the thermodynamic properties on the downstream of the shock wave. 

In the present work, the equilibrium thermodynamic properties for air have 
been obtained from the curve fits of Tannehill and Mugge (ref. 7). 

From the conditions on the downstream side of the shock wave, the fluid 
on the stagnation streamline compresses isentropically until the stagnation 
point is reached. The conditions at the stagnation point can be obtained by 
integrating the following set of equations; 

^ = -pV ^12) 

dV ^ 

h = H - V^/2 (13) 

along with equation (11) from conditions immediately behind the shock wave 
(psw> Psw> ^sw> Vsw) to the point where V = 0 (the stagnation 

point). The final results of this solution will produce the stagnation point 
conditions (pg, Pg, Tg , and hg = Hg). 


Dissociation Enthalpy 

For Lewis numbers other than 1, calculation of the stagnation point 
heating rate from equation (7) requires the dissociation enthalpy for the 
fluid at the edge of the boundary layer at the stagnation point which is given 
by 


(^ 0)3 “ ™(0)(hjj)Q + m(N)(hjj)jj 


(14) 


The dissociation energies ((h|))o and (hj))fj) are readily available, but 
the mass fractions (m(0) and m(N)) must be determined for the stagnation 
conditions of interest. 

The mass fractions of atomic oxygen and atomic nitrogen can be determined 
from the mole fractions as follows 


16 X(0) (15) 

^ ^ 16 X(0) + 32 X(02) + 14 X(N) + 28 X(N2) 

_/xTN _ 14 X(N) (16) 

16 X(0) + 32 X(02) + 14 X(N) + 28 X(N2) 
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Hansen (ref. 8) gives the following approximate equations for determining the 
mole fractions (neglecting ionization) 


x(0) = 2ei/Z 

(17) 

x(02) = (0.2 - ei)/Z 

(18) 

x(N) = 2 e2/Z 

(19) 

x(N2) = (0.8 - S2)/Z 

(20) 

Z = 1 + + G2 

(21) 

^ -0.8 + Vo*64 + 0.81 (1 + 4pg/Kpi) 

(22) 

2(1 + 4Ps/Kpi) 



C2 = 


-0.4 + y/O.16 + 3.84(1 + 4pg/Kp2) 

2(1 + 4ps/Kp2> 


(23) 


Values of £n(l^]^) and £n(Kp 2 ) are given in Table III of reference 8 as 
functions of temperature. 

The value of the dissociation enthalpy obtained by the procedures 
outlined in this section are approximate but should be sufficiently accurate 
for the heating calculations described in this paper. 


RESULTS AND DISCUSSION 


The ratio of the heating rate at a general three-dimensional stagnation 
point ratloed to that at an axisymmetric stagnation point can be determined 
from equation (1) as a function of the parameter k which is the ratio of the 
principal radii of curvature (Rjj/R^). The ratio (qg/^s^AXI 

plotted as a function of k in figure 3. In the discussion that follows, it is 
assumed that is held fixed and k is varied by varying the value of Rg. 

It is easily seen that the heating rate is the highest when the principal 
radii of curvature are equal (i.e., k = 1). Thus the largest value of heating 
occurs at an axisymmetric stagnation point. As R 2 increases, k decreases 
and the stagnation point heating decreases. For k = 0, a two-dimensional 
stagnation point, the heating rate is approximately 70 percent of the value at 
an axisymmetric stagnation point. 
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While no detailed calculations are known to be available for the heating 
at a general three-dimensional stagnation point, it is possible to make 
detailed calculations for the limiting cases of a sphere and a cylinder. 

These calculations have been performed for several points along the entry 
trajectory of the second space shuttle flight (STS-2). The points selected 
are listed in Table I and are the same as those used in the heating analysis 
of reference 9. The calculations were performed for the following set of 
fixed conditions using the transport properties from reference 8. 

T„ = 2000 °R 

Pr = 0.72 

Rx = 1 ft. 

Le = 1 and 1.4 

The heating rates obtained from the present approximate method are compared 
with heating rates obtained from detailed boundary-layer calculations that 
were made by the present author using a method similar to that described in 
reference 10. 

The results for the sphere are tabulated in Table II. It should be noted 
that at the higher Mach numbers the heating rates calculated by the present 
method for a Lewis number of 1.4 are higher than those for a Lewis number of 
1.0. The difference is due approximately to the dlffusional contribution to 
the heating because of the presence of relatively large amounts of dissociated 
oxygen and nitrogen at these conditions. The value of 1.4 was chosen for the 
Lewis number because it has been used extensively by other investigators (for 
example, refs. 2 and 3). At lower Mach numbers where dissociation is less 
Important, the dlffusional contribution to heating decreases. Results for the 
present method for a Lewis number of 1.4 are compared with the boundary layer 
calculations (which Include dlffusional effects) in figure 4. For the portion 
of the STS-2 entry trajectory considered, the present approximate method is in 
good agreement with the boundary layer solutions. Results for the cylinder, 
which are listed in Table III and shown in figure 5, are similar to those for 
the sphere except that the heating rates are lower (as would be expected). 

The present approximate method for calculating the heating rate at a 
general three-dimensional stagnation point has been shown to agree well with 
boundary layer calculations for the special limiting cases of a sphere and 
cylinder. Thus, this method should provide a reasonably accurate method of 
making equilibrium stagnation point heating calculations for "engineering 
type" design and analysis applications. 
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CONCLUDING REMARKS 


An approximate method has been developed for calculating the heating rate 
at general three-dimensional stagnation points during atmospheric entry. 
Although no other experimental or calculated heating rates are known to be 
available for general three-dimensional stagnation points , the present method 
has been shown to agree closely with boundary layer calculations for the 
special limiting cases of a sphere and a cylinder. Thus, it should provide a 
reasonably accurate method of making equilibrium stagnation point heating 
calculations for "engineering type" design and analysis applications. 
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TABLE I. - STS-2 TRAJECTORY POINTS 


CASE 

NO. 

ALT 

KFT 


Voo 

FT/ SEC 

^00 - 

LB/FT^ 

T 

■^00 

°R 

1 

246.0 

25.5 

23620. 

.0454 

357 

2 

219.2 

19.9 

19860. 

.1677 

414 

3 

198.7 

15.7 

16420. 

.3973 

455 

4 

173.8 

11.9 

12680. 

1.0635 

472 

5 

156.4 

9.15 

9710. 

2.0950 

468 

6 

141.5 

7.19 

7540. 

3.8000 

457 

7 

129.8 

6.25 

6460. 

6.1590 

445 
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^APPROX 
BTU/FT^ - SEC 


|btu/ft^ - 
SEC 


Le = 1 


Le = 1.4 


55.7 


63.1 


60.6 


58.9 


65.0 


63.8 


48.4 


52.7 


51.2 


33.8 


36.3 


34.9 


19.0 


19.7 


19.3 


9.87 


9.92 


9.65 


6.35 


6.35 


6.13 























^APPROX 
BTU/FT^ - SEC 

BTU/FT^ - 
SEC 

B 

EBB 

39.4 

44.6 

44.4 

41.7 

46.0 

46.7 

34.2 

37.3 

37.5 

23.8 

25.7 

25.7 

13.4 

14.0 

14.3 

6.98 

7.02 

7.12 





4.44 


4.44 


4.53 






























□o 


O Present method (Le = 1.4) 

□ Boundary layer (includes diffusion) 
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